Thermodynamics of the one-dimensional SU(4) symmetric spin-orbital model 
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The ground state properties and the thermodynamics of the one- dimensional SU(4) symmetric spin 
system with orbital degeneracy are investigated using the quantum Monte Carlo loop algorithm. 
The spin-spin correlation functions exhibit a 4-site periodicity, and their low temperature behavior 
is controlled by two correlation lengths that diverge like the inverse temperature, while the entropy 
is linear in temperature and its slope is consistent with three gapless modes of velocity 7r/2. The 
physical implications of these results are discussed. 



In many transition metal oxides, the electron configu- 
ration has an orbital degeneracy in addition to the spin 
degeneracy. The sign and magnitude of the spin-spin 
interactions is then determined by the orbital occupa- 
tion leading to strong coupling between orbital and spin 
structure (for an overview see Ref. [jjj). The Hamiltonian 
describing such spin-1/2 systems with two-fold orbital de- 
generacy (isospin r = 1/2) was derived by Castellani and 
coworkers more than 20 years ago M . The Hamiltonian 
has rotation symmetry in S'-space. In r-space this sym- 
metry is broken by a Hund's rule term. Recently, the 
investigation of these spin-orbital models has attracted 
renewed interest, following the progress in the experi- 
mental studies of transition metal oxides |3j . 

In this Letter we study the Hamiltonian derived by 
Castellani et al. on a ID chain, but neglecting the Hund's 
rule term. In this isotropic case the Hamiltonian is: 



H = J 



2Si ■ Si+i + — 
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It is rotationally invariant not only in S'-space, but also 
in r-space. Furthermore it has an interchange symmetry 
between spins and orbitals. In such a case, the standard 
mean-field approach || that leads to ferromagnetic cor- 
relations for one type of variables and antifcrromagnetic 
correlations for the other one, should not be appropriate. 
Our main motivation is to study the consequences of this 
symmetry, in more detail. 

A number of analytic results have already been ob- 
tained on this model. The system considered here (Eq. 
(|l|)), belongs to a class of models which is exactly solv- 
able in one dimension by the Bethe Ansatz. The Bethe 
Ansatz solution obtained by Sutherland gives the exact 
ground state energy and the "spin wave" excitations as 
well Q. For the model of Eq. ([!]), there are 3 gapless 
modes, having all a common velocity v — ttJ/2. They 
are shown in Fig. 3 of Ref. 

Second, it was pointed out very recently ||[3) , that the 
Hamiltonian H has not only the obvious SU(2)xSU(2) 
symmetry, but that the full symmetry of Eq. (|lj) is the 



even higher symmetry group SU(4). SU(iV) symmetric 
models in one dimension were studied by Affleck, us- 
ing conformal field theory 0. He showed that any one- 
dimensional system of SU(iV) symmetry is critical. He 
calculated explicitly the critical exponents and zero tem- 
perature correlations and showed that at the very low 
energy scale these models are equivalent to TV — 1 free 
massless bosons. These general results naturally also ap- 
plies to our case with N — 4. 

In this Letter we present the first investigation of the 
thermodynamic properties of the model Eq. (|l|) . For this 
purpose we have adapted the continuous time quantum 
Monte Carlo (QMC) loop algorithm Q to spin-orbital 
models. In fact, this algorithm is so powerful that we 
can also use it to investigate the zero-temperature prop- 
erties of the model: Systems of length L — 100 (with 
periodic boundary conditions (PBC)) and inverse tem- 
peratures (3 J = 400 or 800 (/?>!/);) are predominantly 
in the ground state, and the small contributions from 
thermally excited states negligible compared to our sta- 
tistical errors. In this way, the ground state properties 
can be investigated. We start with a brief summary of 
these results since some of them differ significantly from 
the density matrix renormalization group (DMRG) re- 
sults reported in Ref. ||. 

The ground state energy for a chain of L = 100 with 
periodic boundary condition is found to be eo(L = 100) = 
0.8253(1), in perfect agreement with the Bethe Ansatz 
result for the infinite chain (-0.8251189...) g. The 
zero-temperature correlation functions Wij (T = 0) = 
(SfS?}(T = 0) as a function of \i - j\ (for L = 100) are 
shown in Fig. |l|ab and its Fourier transform in Fig. [l]cd. 

Note that according to the SU(4) symmetry, all the 
following correlations are equal || : 



(S?S?) = (t?t?) = (4S?S?tM) = 



(2) 



independent of the indices a, (3 — x,y, z. This relation 
is valid for zero as well as for finite temperatures. While 
the first equality also holds for an arbitrary SU(2) xSU(2) 
symmetric model with exchange symmetry of the S and 



1 



f-variable, the second one is a special property of the 
SU(4) symmetric model. All the QMC results have been 
checked for the symmetry relation Eq. (^) and perfect 
agreement within the statistical error has been found. 
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FIG. 1. (a) QMC results for the correlation function 
Wij = (SfSj) (solid points) as a function of \i — j\ for a 
chain of length L = 100 with PBC which is predominantly 
in the ground state (for details see text). The correlations 
for \i — j\ = 1,2 and 4 (which are out of the plot range) 
are -0.07168(1), -0.04011(1) and 0.008261(4), respectively. 
Fig. (b) shows the correlations Wij at large distances \i — j\ 
and the fit (for details see text) to the QMC data (open di- 
amonds). The statistical error bars of the QMC calculations 
are much smaller than the symbols, (c) and (d) show the 
Fourier transform S z (k) of Wij on two different scales. 

The correlation function uiij shows a clear four-site pe- 
riodicity (see Fig. [l]). Its sign is positive if \i—j\ = 47V, N 
integer and negative otherwise. The reason for the latter 
is the tendency for every four neighboring sites to form a 
SU(4) singlet ||. Furthermore, from Fig. [jj it can be seen 
that the correlations for distances \i—j\ = 4iV and AN +2 
decay much slower than for \i — j\ = 4N + 1 and 4N + 3. 
The explanation of this fact is simple: The system con- 
sidered here has low lying excitations at k — 0, tt/2 and 
7T (see Fig. 3 of g) each of them leading to a mode with 
wave vector k in the long distance correlations. The am- 
plitudes of this modes are expected all to decay according 
to a power law, but with different critical exponents ctfc. 
From the results for Wij (Fig. [I]), it can be concluded that 
the two dominant modes are those with k = tt/2 (positive 
prefactor) and k = (negative prefactor). This is also 
reflected in the Fourier transform S z (k) of the correla- 
tion function Wy, having a characteristic cusp structure 
at k = 0, tt/2 and tt (see Fig. |l]cd). While the cusps at 
k = and tt/2 are quite sharp, the one at k = n, how- 



ever, is not so pronounced, indicating that the k — tt 
mode is of all the three the least dominant mode in the 
correlation function. 

The two critical exponents 0.^/2 an d «o can be deter- 
mined from the QMC data of the real space correlation 
function w(r) = Wij t u_j\= r . Fitting w(r) to the form 
brr^r-^/z + iL-ry^/^cos^+boir-^ + iL-r)-* ) 
for the range 20 < r < 50 (making explicit use that our 
system has PBC), we find 

a„ /2 = 1.50 ± 0.01, a Q = 1.85 ±0.16. (3) 

The best fit is obtained for b w / 2 = 0.091, a T / 2 = 
1.499, b = -0.035, a = 1.85 and is shown in Fig. 0. A 
precise estimate of ao is not simple since the k = mode 
is only a relative small superposition on the top of the 
much stronger k = tt/2 mode. The exponent a„/2, how- 
ever, can be determined to high precision. These results 
are in very good agreement with the prediction of Affleck, 
who calculated the critical behavior of the SU(4) correla- 
tion function in an arbitrary SU(4) symmetric model us- 
ing conformal field theory 0. This correlation function 
is proportional to Wy, as a consequence of the symme- 



try relation Eq. 



and the exact results are a 



tt/2 



and ao = 2. The exponent awa has also been estimated, 
using DMRG (a ff/2 ~ 1.5 - 2) [§. The DMRG results 
are in principle more precise than the QMC results, but 
finite size effects in DMRG studies are much bigger due 
to the use of open boundary conditions. Thus it is not 
surprising that our estimate Eq. (||) is much more precise. 

At finite temperatures, the dominant components in 
the correlation function, Wij(T) = (SfSj)(T) (note that 
Eq. (||) holds also at finite T) which result from the 
soft modes at k = and 7r/2, no longer decay according 
to a power law, but exponentially. The corresponding 
correlation lengths £o{T) and ^/^{T) may be different. 

The correlation function (SfS^)(T) is shown as a func- 
tion of \i — j\ in Fig. |^ for a system of length L = 200 
with PBC at a temperature T = 0.05 J. To find the cor- 
rect low-temperature form, describing the long distance 
behavior (\i — j\ 3> ^0,^/2) of the correlations Wij(T), 
one has to consider not only a phase shift S(T) in the 
k = tt/2 mode, but also an incommensuration effect of 
this component, i.e. that the period is shifted away from 
k = tt/2 by an amount 4>k(T). This is due to the asymme- 
try of the excitation spectrum at the point k = tt/2 which 
manifests itself at finite T, where also excited states con- 
tribute to Wij(T). This asymmetry can be seen in Fig. 2 
of |(|, where the degeneracy of the lowest "spin wave" 
branch is indicated. As the degeneracy for k > tt/2 is 
larger than for k < tt/2, we expect the weight of the 7r/2 
mode to be shifted to a higher k value. This effect can 
also be observed in the Fourier transform S z (k, T) of the 
correlation function, where the maximum at k = tt/2 at 
T = moves to higher k- values when T increases. 

Finally, we propose the following low temperature form 
for the correlations Wij(T) with \i — j\ 3> Co,Ctt/2 : 



2 



Wij (T) = b (T)e~^^ + 6 x/2 (T)e-l i - i l/^/^ T ) 

cos[(7r/2 + fc (r))(i-i) + *(T)]. (4) 

Fitting the above form to the QMC data of itty (T) for 
various T gives the temperature dependence of the six 
parameters bo, £oj fn-/ 2 ' 0fc and S. The correspond- 
ing fit for T = 0.05 J is shown in Fig. ||. At this point, 
we want to emphasize that it is important to include the 
effect of incommcnsuration (i.e. including a parameter 
4>k 7^ 0) to get accurate fits in the temperature range 
0.01J < T < 0.08J and that at finite temperatures and 
larger distances, also the correlations at distances 4iV + 3 
can become positive (Fig. ||), different from the corre- 
lations at T = 0. For all considered T, the correlation 
lengths are much smaller than the length of the system 
(L > 6^/2, 6£o)j so that finite size effects are negligible. 
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FIG. 2. QMC results for w^{T) = (S|5|)(T) for a system 
of length L = 200 with PBC and at temperature T = 0.05J 
(solid points). As a guide of the eye, the correlations at dis- 
tances \i — j\ = AN + m are connected by a separate (dotted) 
line for each m = 0, 1, 2 and 3. The inset shows the long dis- 
tance correlations. The error bars are much smaller than the 
symbols. The open squares show the fit of form Eq. to 
the QMC data. 

In Fig. ||the inverse correlation lengths t;^ 1 and £~L as 
well as the periodicity shift are plotted as a function 
of temperature. At very low temperatures, ^q 1 and £~J 2 
both show a linear behavior (see Fig. ||) and the leading 
temperature dependence are found to be 

= (2.99 ± 0.03)T, f- 1 = (3.90 ± 0.09)T (5) 

Therefore both correlation lengths scale with 1/T. This 
scaling behavior, including the prefactor, can be moti- 
vated in the following way. By the Lorentz invariance of 
the underlying field theory of the considered model and 
by exchange of the imaginary time and space direction, 



one has ( k (T) = v/A k {L = v/T), where A k (L) is the fi- 
nite size gap to the lowest excitations at wave vectors « k 
in a system of length L. v is the spinon velocity, which in 
our model is irJ/2. For k ss 0, the lowest lying excitation 
energy is A (L) = v ■ 2tt/L leading to £o = J/(4T), in 
good agreement with Eq. (|J)). For k = tt/2, the finite 
size results of Ref. |6) show that A^/ 2 {L) ~ 0.75A (L) 
(a similar result should be obtained, using field theory) 
and hence ~ J/(3T), again in very good agreement 
with Eq. (Q). In comparison, the leading temperature 
dependence of the correlation length in the SU(2) AF 
Heisenberg model (iJ H B = JY^i^i' w hb = ttJ/2 

and A™(L) = irv/L) is £™(T) = J/{2T). 

The leading T dependence of the periodicity shift 4>k is 
determined by fitting XT 13 to the data at very low temper- 
ature <pk{T) oc J 1 2-ii±0.i5_ Thjg scaling exponent is quite 
close to the value of 2, which one would expect from a 
simple calculation, considering the thermal admixtures 
of the spin wave branches. 




FIG. 3. Inverse correlation lengths £ 1 (open circles) and 
Cr/2 ( n ll e< i circles) and the periodicity shift (f>k (inset) as a 
function of temperature. The leading T dependences (Eq. 
(^)) are also shown (solid, dotted and dashed lines). 

Finally, we concentrate on the entropy s per site of 
the SU(4) invariant model Eq. (|l|). Its T-dependence 
s(T) is shown in Fig. |^. With decreasing T, the en- 
tropy decreases monotonically from the high temperature 
value In 4 to at zero temperature. At low temperatures 
the entropy shows a linear behavior as in the AF SU(2) 
Heisenberg chain (Hub ) ■ The slope in the spin orbital 
model, however, is about a factor three bigger than that 
in the AF Heisenberg chain (see inset of Fig. ^). This 
is consistent with the statement of Affleck that the 
AF Heisenberg model is equivalent to one free massless 
boson, while the SU(4) invariant spin orbital model is 
equivalent to three massless bosons. The velocity of these 
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bosons are all equal to irJ/2 [Q,|9|. Therefore we expect 
the low energy density of states (and hence the entropy) 
of these two models just to differ by a factor 3. 

The implications of these results for mean field treat- 
ments are far reaching. To put them in perspective, it is 
useful to compare them to the standard mean-field de- 
coupling Q (Si ■ Si+i)(fi ■ f i+ i) -* < Si ■ S i+ i > Ti ■ f i+ i 
+ < n-n+i > S l -S l+ i - < Si-S i+ i X n-n+i >. Such 
a decoupling has a number of consequences. First of all, 
the correlation function ((Si ■ Si+i)(fi ■ 7*4-1)) should be 
equal to the product of < Si ■ Dj+i > with < 7* ■ 7*4-1 >, 
in clear contradiction both with the fact that all of them 
are negative according to our results and with the prop- 
erty of Eq. (3). Besides, and more importantly, if such a 
decoupling was a valid approximation, the low-lying ex- 
citations should consist of two branches corresponding to 
spin and orbital excitations, respectively. This is again 
in clear contradiction with the 3 low-lying modes of the 
Bethe ansatz which control the low temperature physics 
according to our entropy results. So there is a manifest 
break-down of the mean-field decoupling when spin and 
orbital degrees of freedom play a symmetric role. 
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FIG. 4. Temperature dependence of the entropy s per site 
for the spin-orbital model Eq. (|l]) (solid line). In the inset 
the entropy per site is shown on larger temperature scale to- 
gether with the entropy shb per site of a SU(2) spin-1/2 AF 
Heisenberg chain (Hub) (dotted line). For comparison also 
3shb is shown (dashed line). 

What is then the nature of the low-lying excitations? 
A full answer cannot be given on the basis of the present 
results, but a number of conclusions can be reached. Let 
us start with the 2-site problem. The ground state is 
six- fold degenerate (spin-triplet x orbital-singlet or spin- 
singlet x orbital-triplet), and energy may be gained by al- 
lowing fluctuations between these local configurations. 
The mean-field decoupling fails because it cannot take 
advantage of these fluctuations. Elementary considera- 



tions show that the best mean-field decoupling leads to 
a very poor estimate of the ground state energy (-0.3863 
vs -0.8251 for the exact result). That it is possible to 
gain energy by allowing the system to fluctuate locally is 
best exemplified by the 4-site problem. In fact, the ex- 
act ground state for a 4-site cluster with periodic or open 
boundary conditions, the SU(4) sing let ofRef. @, can be 
obtained exactly in terms of these dimer wave functions, 
and the energy per bond is equal to -1, i.e. each bond 
has now managed to reach its ground state energy thanks 
to the fluctuations between these 6 local configurations. 
Note that this is no longer true for longer systems, indi- 
cating that 4-site clusters should be a good starting point 
for building variational wave- functions. This can be seen 
as the physical origin of the 4-site periodicity of the cor- 
relation functions. A similar conclusion was reached in 
Ref. H on the basis of the SU(4) symmetry. 

Finally, if the ground state is an RVB-like state involv- 
ing resonances between different local configurations, we 
are lead to the conclusion that the elementary excitations 
cannot be pure spin or orbital excitations, but compos- 
ite objects where spin and orbital degrees of freedom are 
intimately connected. Work is in progress to get a more 
precise picture of these excitations. In addition, the pre- 
sented results will also have dramatic consequences for 
more realistic models where the interchange symmetry 
between spin and orbital degrees of freedom is only ap- 
proximately valid. This is left for future investigation. 
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